Screening effects in superconductors 
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The partition function of the Hubbard model with local attraction and long range Coulomb repulsion 
between electrons is written as a functional integral with an action A involving a pairing field A and a 
local potential V . After integration over V and over fluctuations in |A| 2 , the final form of A involves 
a Josephson coupling between the local phases of A and a "kinetic energy" term, representing the 
screened Coulomb interaction between charge fluctuations. The competition between Josephson 
' coupling and charging energy allows to understand the relation between Tc and composition in 

^\ , high- Tc materials, in particular superlattices, alloys and bulk systems of low doping. 
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We start from a Hamiltonian describing charge carriers (electrons or holes) on a lattice, subject to an on-site 
■ attraction, — U, and the long range Coulomb repulsion, e 2 Vc, acting on particles on different sites : 
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Here (c ; a ) are creation (annihilation) operators for charge carriers with spin a at lattice site I, {I, I') denotes pairs 
O ■ of nearest neighbor sites, \i is the chemical potential, uq the background, neutralizing the density of charge carriers, 
and Ni = c / + <r c / o-- The partition function can be written as a functional integral by means of two successive 
Stratonovich-Hubbard transformations decoupling the two interaction terms in H with the help of a complex 

field A and a real field V : 
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We then evaluate the trace over the electronic degrees of freedom 



and we expand the "free energy" F up to fourth order in A and to second order in V and up to leading terms in 
space and time gradients of the two fields : 
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Here the coefficients a, b, c, d are related to the free electron particle-particle propagator is the free electron 

contribution, and Vg C (l, I') = Vq (I, l')+Xo(l, I') is the screened Coulomb potential (approximated by its static limit) 
with Xo(M') being the electronic polarizability. Integrating over the electric potential V yields 
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where p(l,r) — 2d\A(l,r)\ 2 + (Ni) — no represents charge density fluctuations. Next we introduce amplitude and 
phase of A. In the following, we assume to be in a relatively strong coupling regime (t <C U) in which fermions bind 
into onsite singlet pairs at a temperature on the order of the mean field transition temperature T^j, well above the 
superconducting phase transition, the latter being finally triggered by the onset of phase order Below T m f, 

a < 0, so that the average amplitude has a non- vanishing mean value Aq given by a + 2b A 2 ,. Charge neutrality implies 



2dA 2 + (Ni) -n Q = and, for t < U @|-|§ 



A 2 = no (2- n ) Tj2 _ n Jj2 
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for n < 1. 
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(neglecting gradient terms) yields a free energy functional for phase fluctuations only 



Splitting the number of pairs at a given site into |A(Z,r)| 2 = A 2 , + ^-Sn p (l,T), we integrate over Sn p (l,T) which 
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with the Josephson coupling [jiJU J — 2cAq = 2n ( , and 
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Our expansion of F in powers of A has yielded the on-site repulsion (2e) 2 ^ between pairs. However, due to the 
exclusion principle, two pairs cannot really sit on the same lattice site. Thus, in the following, we exclude I = I' in 
the last term of (13). By going from the "phase velocities" ^ to the conjugate momenta p(l, r) we end up with the 
partition function of the Hamiltonian ]10[ | 



H = \ E (K0 - y) [i^) 2 w(i, 0] _ !£) + j £ [i - cos (6(0 - e(0)]. 
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Hamiltonian (15) describes the relevant physics of short coherence length superconductors in terms of Josephson 

coupled spatial phase variations and " charge fluctuations" coupled by the screened Coulomb interaction. It is also the 

" phase-only" representation of the Hamiltonian of interacting bosons Jl0|-[r2[ . Its critical behavior, in particular the 

influence of the "background" , — has been studied in these references. Here, we apply expression (15) to calculating 

the transition temperature of strongly anisotropic superconductors, such as superlattices and bulk systems in the 

underdoped regime OL We make the following approximations : i) H is restricted to one superconducting layer; 

ii) the screened Coulomb interaction, which takes into account the electric coupling between layers, is modelled by a 

Yukawa-form, with a Thomas- Fermi screening length Atf, depending on the density no of charge carriers according 

i 

to the wellknown formula Atf = ^ [A -1 ], where r s = ^ 4 ^ ng J 3 and ao is the Bohr radius; iii) considering 

only no <C 1, the "background shift", — in the first term of (15) is neglected. We make connection with previous 
work fl4|,|l5| by mapping (15) onto a "capacity model": 

i$>(0 [{2efW{l,l')]p{l') « lg!j>(0 2 , (16) 
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with ^ = ^J2i W(l, 0) = e -^F . 

This is the XY-model with kinetic energy |l4|,[l5|], e being the dielectric constant of the interlayer material and 
the lattice constant. In Ref. [14] the critical temperature in two dimensions was evaluated in the "self-consistent 



harmonic approximation" (SCHA) which gives results in good agreement with Monte Carlo simulations |16|. A good 



overall fit of the numerical SCHA result is 17 Tc{a) w ?c(0). /l — where a = ^L. is the ratio between charging 



a c ""^^ " — 2CJ 

and Josephson energy. When a approches ac=6.2, Tq goes to zero. This approach has been successfully applied 
|u],(l5| to calculating the superconducting transition temperature for superlattices and for alloys by determining the 
appropriate effective capacity through electrostatic considerations. 

We finally use Hamiltonian (15) to find Tc in function of doping for high-Tc superconductors in the underdoped 
regime. For YBa2Cu307, our "capacity model" has allowed to fit the Tc variation of both, superlattices and alloys, 
in a coherent way, using SCHA Jl4|,[l5| for J—12UK at optimal doping (no ~ 0.16 of holes per cell (ljj). Below optimal 
doping, the ratio a increases when the number no of charge carriers is reduced : the no-dependence of J is given 
following Eqn. (13) and C varies with no through the screening length Atf- Using the above expression for Tc{a), 
with doping-dependent J and C, we then find that Tc should be zero for no ~ 0.07, in good agreement with the 
measured phase diagram This shows that the phase boundary in the underdoped regime can be understood 

in terms of Bose-Einstein condensation of preformed pairs, Tc being suppressed by phase fluctuations, when the 
minimum doping is approached. 

In summary, starting from the attractive Hubbard model with long range Coulomb repulsion, we have given a micro- 
scopic derivation of a description of short coherence length superconductors in terms of the superconducting phase, 
the Hamiltonian for which includes a "charging energy" and a Josephson coupling. This is a microscopic justification 
of such a Hamiltonian, which has been used previously |l4|,|l5| to calculating Tc for superlattices and alloys and also 
allows understanding the phase boundary of bulk oxides in the underdoped regime. 

This work was supported by the Swiss National Science Foundation. We thank T. Schneider and M.H.Pedersen for 
valuable discussions. 



[1] T.K. Kopec and H. Umezawa, Phys. Rev. B 47, 8923 (1993). 

[2] V. Ambegaokar, U. Eckern and G. Schon, Phys. Rev. Lett. 46, 211 (1982). 

[3] M. Drechsler and W. Zwerger, Ann. Phys. 1, 15 (1992). 

[4] M. Randeria, in Bose Einstein Condensation, ed. by A. Griffin, D. Snoke and S.S. Tringari, Cambridge University Press 
(1994). 

[5] M.H. Pedersen and H. Beck, not published. 

[6] V.J. Emery and S.A. Kivelson, Nature 374, 434 (1995) and Phys. Rev. Lett. 74, 3253 (1995). 

[7] S. Doniach and M. Inui, Phys Rev. B 41, 6668 (1990). 

[8] E. Roddick and D. Stroud, Phys. Rev. Lett. 74, 1430 (1995). 

[9] J.J. Vicente Alvarez and C.A. Balseiro, Solid State Comm. 98, 313 (1996). 
[10] MP. A. Fisher and G. Grinstein, Phys. Rev. Lett. 60, 208 (1988). 
[11] Min-Chul Cha and S.M. Girvin., Phys. Rev. B 49, 9794 (1994). 

[12] M.P.A. Fisher, P.B. Weichmann, G. Grinstein and D.S. Fisher, Phys. Rev. B 40, 546 (1989) and A. van Otterlo, K.H. 
Wagenblast, R. Baltin, C. Bruder, R. Fazio and G. Schon, Phys. Rev. B 52, 16176 (1995). 



3 



[13] J.L. Tallon and N.E. Flower, Physica C 204, 237 (1993) 
[14] D. Ariosa and H. Beck, Phys. Rev. B 43, 344 (1991) 

[15] D. Ariosa, T. Luthy, V. Tsaneva, B. Jeanneret, H. Beck and P. Martinoli, Physica B 194-196, 2371 (1994). 
[16] L. Jacobs, J.V. Jose, M.A. Novotny and A.M. Goldman, Phys. Rev. B 38, 4562 (1988). 
[17] D. Ariosa, J. Perret, H. Beck and P. Martinoli, in preparation. 



4 



